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■ ABSTRACT 
O " 

o 

CN ' We present Keck LWS images of the Orion BN/KL star forming region ob- 

• tained in the first multi- wavelength study to have 0('3-0!'5 resolution from 4.7 fim 

[T I ■ to 22/im. The young stellar objects designated infrared source n and radio source / 

00 ■ are believed to dominate the BN/KL region. We have detected extended emission 

from a probable accretion disk around source n but infer a stellar luminosity on 
the order of only 2000 Lq. Although source / is believed to be more luminous, we 
do not detect an infrared counterpart even at the longest wavelengths. However, 



> 



CN ■ we resolve the closeby infrared source, IRc2, into an arc of knots ~ 10^ AU long 

^ ■ at all wavelengths. Although the physical relation of source / to IRc2 remains 

Q . ambiguous, we suggest these sources mark a high density core (10^-10^ pc~^ over 

O ■ ~ 10^ AU) within the larger BN/KL star forming cluster. The high density may 

^ . be a consequence of the core being young and heavily embedded. We suggest the 



o 



energetics of the BN/KL region may be dominated by this cluster core rather 
than one or two individual sources. 

Subject headings: infrared: stars — ISM: individual (Orion KL) — stars: forma- 
tion — stars: pre-main-sequence 



1. Introduction 

Orion BN/KL is the closest region of high-mass star formation (~ 450 pc) and though 
also the archetype, it is not necessarily well understood (e.g., Genzel & Stuzki 1989). For 
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instance, the sources that drive extensive outflow activity in BN/KL (e.g., Wright et al. 1996) 

and contribute most to its luminosity have not been identified with certainty (e.g., Menten & 
Reid 1995). This is true because BN/KL is crowded, exhibiting ~ 20 compact infrared (IR) 
peaks distributed over ~ 20" and because extinction is high and variable across the region 
(Gezari et al. 1998). As high-mass stars often form in clusters (Garay & Lizano 1999, and 
references therein), there are presumably numerous embedded protostars and young stellar 
objects (YSOs) in the region, many of which have probably not been recognized yet. 

Two notable sources in BN/KL exhibit centimctcr-wave continuum and as a result 
are believed to be luminous YSOs, IR source n (Lonsdale et al. 1982) and radio source / 
(Churchwell et al. 1987). Source n, also known as radio source L (Churchwell et al. 1987), 
exhibits a bipolar radio jet (Menten & Reid 1995), hard X-ray emission (Garmire et al. 2000; 
Feigelson et al. 2002), and A2.3//m CO overtone emission probably from the hot surface of 
an accretion disk (Luhman et al. 2000). Source /has been difficult to study because it has 
not been detected at IR wavelengths, probably due to high extinction. It lies close to the 
mid-IR source IRc2 (Gezari 1992; Menten & Reid 1995) for which Ay ~ 60 (Gezari et 
al. 1998) and the peak of the Orion Hot Core, for which probably Ay ^ 10^ (Genzel & 
Stuzki 1989; Migenes et al. 1989; Plambeck et al. 1995). However, source /excites II2O and 
SiO maser emission in dense, hot material at radii < 10'^ AU (Genzel & Stuzki 1989; Gezari 
1992), which may be used to study gas dynamics in detail (e.g., Greenhill et al. 1998, 2003). 

We have imaged BN/KL in the mid-IR at resolutions of 0''3-0"5 with the intention to 
characterize sources of excitation. This is the first multiwavelength study with sub-arcsecond 
resolution up to 22/im. (cf., Dougados et al. 1993; Rouan et al. 2004; Shuping et al. 2004). 
We concentrate on a largely qualitative interpretation of the data. Quantitative modeling 
and analysis of the whole region will be discussed in a separate paper. 

2. Observations and Reduction 

We observed the BN/KL region with the Keck-I telescope on 1999 December 21, 2000 
October 14, and 2002 August 20. We used the Long Wavelength Spectrometer (LWS; Jones 
& Puetter 1993) in imaging mode at seven filter wavelengths from 4.7//m to 22.0/im (Table 1) 
to obtain frames of 128 x 128 pixels (~ 10"x 10"), with which we created mosaics. We chopped 
at 5 Hz and nodded the telescope both using fields 30" west in 1999 and 2002 and 40" north 
in 2000. Registration of overlapping frames is uncertain to < 1 pixel. To register images 
obtained at different wavelengths, we used the peak of the Becklin-Neugebauer object (BN) 
as a reference. When BN lay outside a mosaic, we used IRc7. The formal uncertainty in 
registration is < 1 pixel. 
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We used observing procedures and data reduction techniques discussed by Gezari et 
al. (1992) and implemented in the MOSAIC image processing package (Varosi & Gezari 
1993). Observations of a CMa demonstrated the high quahty of the point source responses. 
Images at 12.5/im and longer wavelengths are diffraction limited (0''3-0''5), while seeing, 
chopping/nodding error, and smoothing during processing limit resolution to ~ 0''3 for 
shorter wavelengths. To calibrate flux densities, we scaled images at each wavelength to 
match the intensities for BN or IRc7 reported by Gezari et al. (1998), interpolated to the 
LWS passbands (Table 1). 

We calibrated the pixel scale to first order using observations in 1999 of the binary 
XYPer for which we adopted a 1"33 separation and 256.3° position angle (Bertout, Ro- 
bichon, & Arenou 1999). We made second order adjustments to align the IR and radio 
centroids of BN and source n on a A12.5//m mosaic, where we defined the radio centroid for 
source n as the midpoint between the lobes. The adopted scale and rotation were 0''0840 
pixel"^ and 2.1° clockwise (i.e., north is rotated counterclockwise from vertical on the im- 
ages). Uncertainties of 1-2% and 0.5° include the ~ 0''02yr^^ proper motion for BN (Plam- 
beck et al. 1995). We rechecked the plate scale in 2002 with observations of a Her, for 
which we adopted a 4''84 separation and 105° position angle (Mason et al. 2000; Docobo et 
al. 2000), and obtained 0''0831 pixel"^ and 4° clockwise, with uncertainties of 2% and 1°, 
in good agreement with earlier results. 

We adopted the 1999 plate scale for images obtained in 2000 because a parallelogram 
distortion precluded our estimating a new solution. However, because the sub-fields shown 
here are small (containing just IRc2 and IRc7) the resulting errors are < 2 pixels. We 
estimated this error empirically by registering the peak of IRc7 observed at A12.5//m in 2000 
and 1999, applying the registration offsets to images from 2000 for other wavelengths, and 
measuring the misalignment in the peaks of IRc2. 

To register the radio position of source / on the IR images, we assumed the radio and 
IR centroids for BN were coincident and used the angular separation of BN and source/ 
estimated by Menten & Reid (1995). For images obtained in 2000, we measured angular 
offsets from (closeby) IRc7 to minimize the effect of the parallelogram distrotion. The total 
uncertainty in radio-IR registration is 1-2 pixels. (We note that the proper motion of BN 
contributes only a ~ 0.4 pixel error to the registration.) 
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3. Results and Discussion 

3.1. Source n as an intermediate luminosity star 

We have found that emission from source n at XS.O/im and A11.7//m is elongated nearly 
east-west (Figure 1) and orthogonal to the bipolar radio structure detected previously by 
Menten & Reid (1995). In contrast, other sources in the A8.0/xm and A11.7//m images are 
not elongated (e.g., lRc7). We propose that the elongation arises from warm dust in the 
accretion disk inferred previously from the presence of CO overtone emission at \2.3fim 
(Luhman et al. 2000). We fit a single component elliptical Gaussian model to each image. 
At A8.0;um, the source size (half power full width)is (0''74 x O^'SO) ± 0''08, and at A11.7/im, 
the size is (0^75 x 0'.'53) ± 0:'12 (~ 340 x 230 AU). The position angles are 115° ± 13° and 
101° ± 24°, respectively, while the position angle of the bipolar radio structure is 12 ± 1°. 
The shape of source n at longer wavelengths is more uncertain, because of plateau emission 
between IRc4 and IRc7. However, if we assume the disk emits as a blackbody peaking close to 
the shortest wavelength at which it is detected (8.0/xm), then the characteristic temperature 
is on the order of 300 K. 

Although source n is extended at mid-IR wavelengths, it is point-like in high angular 
resolution observations from A2.2//m to A5//m (Dougados et al. 1993; Stolovy et al. 1998). 
This near-lR emission may arise from the inner edge of the accretion disk, close to where 
dust is sublimated. For the disks of Herbig AeBe stars that are resolved with interferometry, 
this inner radius is < 0.2 - 10 AU (Tuthill et al. 2001; Monnier & Millan-Gabet 2002). 
At the distance of BN/KL, the near-lR emission would subtend < 0''02, thus explaining its 
point-like appearance in observations so far. 

The luminosity of source n is not known. However, mid-lR emission from an accretion 
disk may be expected to arise from surface material heated by stellar radiation, and we can 
use the characteristic radius, R, and temperature, T, to estimate the luminosity of source n, 
L*. Dust near the surface of the disk is optically thick to stellar radiation (mostly optical 
and ultraviolet) but optically thin to re-emission (mostly mid-IR). For thermal equilibrium, 
~ IGTrR^aT^lQe/Qa), where Qe and Qa are grain emissivities for emission and absorption, 
and (7 is the Stefan-Boltzmann constant. For absorption of stellar radiation, Qa — 1, while 
for re-emission, Qe ^ 2'ria/X, where a is a characteristic grain size. For A ~ lO/im and 
a ~ O.l/im (typical for interstellar dust), Qe/Qa ~ 0.06 and ~ 2000 L© (mid B-type). 
This luminosity is consistent with the observation that at A1.6yum source n is ~ 10 times 
dimmer than BN (1-2x10^Lq; Scoville et al. 1983), given that source n is 1.2 magnitudes 
less red between Xl.l/im and 1.6/xm (Luhman et al. 2000), if we assume the two stars 
have approximately the same intrinsic near-IR color. Nonetheless, the estimated luminosity 
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should be regarded with caution. It depends on uncertain values for a and T rather than 
radiative transfer modeling of the spectral energy distribution, and the estimated luminosity 
exceeds the apparent luminosity of ~ 10 (Table 1). Moreoever, identification of source n 
as a hard X-ray source is a concern since a B-type star is unlikely to be active, though the 
X-rays may be associated with a less massive companion or the jet. 

3.2. Source /, IRc7, and IRc2 

An IR counterpart to radio source / is long soughtafter. The LWS observations are the 
first to resolve the angular separation between source / and the comparatively bright IRc2 
at wavelengths up to 22fim. Nonetheless, we do not detect emission from a counterpart 
with an upper limit on the order of 1 Jy at \8.0fim and 10 Jy at A22/im. Wc suggest that 
emission from source /is blocked by the Orion Hot Core. Plambeck et al. (1995) estimate 
the average optical depth of the core is ~ 0.05 at A3. 5 mm. The wavelength scaling of opacity 
to shorter wavelengths is uncertain because it depends on grain composition and the amount 
of processing. However, if we assume spherical silicate dust grains (Mathis 1990), then the 
optical depth is > 300 at A8.0/xm and A22//m. This is sufficient to obscure the mid-IR dust 
emission from a 10^ K blackbody of diameter up to 80 AU and luminosity up to ~ 6 x 10^ Lq, 
which is comparable to the luminosity of the BN/KL region (Gezari et al. 1998). It is also 
sufficiently large to explain nondctcction at \2.2fim in Hubble Space Telecope (HST) and 
Keck observations (Stolovy et al. 1998; Hillenbrand & Carpenter 2000). 

IRc7 is relatively compact and prominent at all observed wavelengths longer than 8/xm 
(Figures 2 and 3). This is probably indicative of star formation at a relatively early phase. 
IRc7 also lies close (in projection) to the radio axis of source n (Stolovy et al. 1998), and 
HST observations at A2.2/im have detected a fan of emission opening away from source n 
with a bright point at the vertex. This morphology and orientation may signify interaction 
with radiation and outflow incident from source n. Alternatively the near-IR fan may be a 
scatting cone along the rotation axis of a protostellar system embedded in an optically thin 
envelope or halo. 

IRc2 is shaped hke a ~ 10^ AU long arc that opens to the southwest, in the direction of 
source/, which hes close to one end (Figure 2). The locations of knots along the arc depend 
on wavelength and are distributed 0''6 to 2"7 from source/. Similar emission distributions 
at A3. 8, 5, and 12.5;um are reported by Dougados et al. (1993), Rouan et al. (2004), and 
Shuping et al. (2004). Emission at the shorter wavelengths is the closest to source /while 
emission at A22.0;um is among the farthest (Figure 3). Although suggestive of central heating, 
emission at wavelengths near the peaks of the silicate absorption features at 9.7/xm and 
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18.5//m also tends to lie away from source / (Figures 2 and 3). If local extinction toward 
source / is enhanced, then it is surprising that the short wavelength emission arises closeby 
to source/. At the same time, we note that the shortest wavelength emission yet detected 
from IRc2, at A2.2/im, is associated with the knot "B" (Dougados et al. 1993), ~ 1" from 
source / (Stolovy et al. 1998; Hillenbrand & Carpenter 2000). 

The relationship of IRc2 to source / is ambiguous. Because massive stars form in clus- 
ters, we speculate that the ~ 5 knots in lRc2 may be individual sites of star formation, which 
together with source / form at least part of a high density core within the larger (~ 20") 
BN/KL star forming cluster. In support of this we note that (1) the relative distributions of 
emission in IRc2 from A2.2/im to A22.0/im do not exhibit an overall systematic pattern (e.g., 
central heating), and (2) knot "C" (Dougados et al. 1993) is both closely associated with a 
cluster of OH masers (Figure 2) and X-ray emission (E. Feigelson, private communication). 
(If the X-ray emission is the result of a wind driven by source/, then other knots would 
probably be X-ray sources as well.) The inferred characteristic density of lO'^-lO^pc"^ over 
~ 10'^ AU is unusually large and greatly exceeds the density of the closeby Trapezium (Hil- 
lenbrand & Hartmann 1998). Because high cluster densities are associated with high stellar 
masses (Tcsti et al. 1999, and references therein) we might expect to observe precursors to 
0-type stars, but the large density may be due to cluster youth. As clusters evolve they 
expand because of stellar dynamics and the expulsion of gas (Kroupa 2004, and references 
therein) . Clusters that are still deeply embedded may be expected to exhibit unusually high 
central densities though observational evidence is weak. 

We cannot yet exclude the possibility that lRc2 is heated externally, by source/. The 
tip of IRc2 recedes systematically as wavelength increases from A5//m to A22.0/xm (Figure 2), 
which is suggestive of a local temperature gradient. In principle, IRc2 could mark the edge 
of a wind-blown cavity. Relying on maps of SiO masers Greenhill et al. (1998) propose that 
source / drives a conical outflow to the north and west, which IRc2 just happens to subtend. 
Greenhill et al. (2003) present an alternative model in which source /binds an edge-on flared 
disk whose limbs extend north and west from the star. In this case lRc2 straddles the disk 
plane. The mid-lR arc would comprise irradiated disk surface material and peripheral low 
latitude material heated by other YSOs. 



4. Conclusions 

We have detected the probable signature of an accretion disk around source n at mid-IR 
wavelengths, and we estimate the stellar luminosity is ~ 2000 L©. We have not identified 
any counterpart to radio source /, despite higher angular resolution and broader wavelength 
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coverage than previous studies. We resolve IRc2 into knots at up to A22.0//m and propose 
that they are sites of star formation. Together, source / and IRc2 may mark the dense core 
(10''-10^ Mq pc^^) of the larger BN/KL infrared cluster. Many past studies have sought 
to demonstrate that luminosity and outflow associated with BN/KL are driven by one or 
two sources. This seems increasingly unlikely. We note that source n has only a modest 
luminosity and its jet axis is poorly aligned with outflow axes observed on scales of 10^ AU 
(e.g., Schultz et al. 1999). Moreover, members of the putative IRc2-source / subcluster are 
ready candidates to drive outflow collectively. In a later paper, we will present a radiative 
model of the LWS data, with estimates of dust temperatures, hne-of-sight optical depths, 
and intrinsic luminosities for the dominant mid-IR sources in BN/KL. 

We thank E. Feigelson, S. Wolk, and the COUP consortium for sharing X-ray source 
coordinates before publication. F. Varosi provided technical support for the MOSAIC soft- 
ware. The Keck Observatory was made possible by the flnancial support of the W. M. Keck 
Foundation and is operated as a scientiflc partnership among CalTech, the University of 
California, and NASA. 
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Table 1. Peak Specific Intensities 



Filter 


Band Ep.(«)lRc2 


IRc7 


n 


BN 


(/xm) 


(/xm) 




(Jy arcsec"^^ 






4.7 (M) 


4.4-5.0 2 


0.64 


0.15 




58 


8.0 


7.5-8.2 1 


15 


5.0 


4.4 




10.7 


10.0-11.4 3 


1.3 


2.6 


0.96 


93 


11.7 


11.2-12.2 3 


4.9 


7.2 


2.4 


140 


12.5 


12.0-13.0 1 


13 


17 


3.3 


152 


17.65 


17.3-18.2 3 


16 


34 


2.3 


175 


22.0 


21.0-23.0 3 


43 


86 


3.2 


190 



("^ Epoch for tlie tabulated photometry (see Sec- 
tion 2). 

'^''^ Intensities are uncertain to ~ 30%. The zero level 
was taken from a sky position ~ 3" south of source /. 
Plateau emission between IRc4 and IRc7 has been sub- 
tracted in the case of source n. Images are scaled to 
match listed IRTF peak intensities for BN, or IRc7 at 
A8.0/im (Gezari et al. 1998) interpolated to LWS pass- 
bands. 




Fig. 1. — Images of source n obtained in 1999. The elongations at A8.0 and A11.7/im are 
within ~ 10° of perpendicular to the jet axis (red bar). The peak at each wavelength is 
listed in Table 1, and the specific intensity scale is linear. The upper right bar in each panel 
indicates north. 



Fig. 2. — Mid-IR images obtained in 1999 and 2000 (A4.7/im only), (right) Mosaic of 
A12.5/im emission. Contours correspond to 1.2, 1.6, 2.1, 2.9, 4.0, 5.6, 7.8, 11, and 16 Jy 
arcsec^^. The white bar indicates north. The box outhnes the region shown in close-up. (left) 
Close-ups of IRc2 and IRc7. Black crosses mark the approximate positions of A12.5;um peaks 
and are provided for registration. Symbol size indicates the alignment uncertainty of the 
panels. The red cross marks the position of source /, and symbol size reflects uncertainty. The 
distributions of H2O (red circles) and OH masers (white squares) are also shown (Greenhill et 
al. 1998; Johnston et al. 1989). The H2O masers lie in a high-density outflow from source/ 
(arrows). The OH masers lie in lower density gas and are probably pumped by far-IR 
radiation. The peak speciflc intensities of IRc2 and IRc7 at each wavelength are listed in 
Table 1, and the speciflc intensity scale is linear. 



- 13 - 




Fig. 3. — Images of IRc2, IRc7, and source n obtained in 2002. Peak specific intensities 
are listed in Table 1, and the specific intensity scale is linear. The white bar in each panel 
indicates north. 



